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Abstract 



One of the most intriguing unsolved problems in modern cosmology involves the 
nature of the observed dark matter in the universe. Various hypothetical particles have 
been postulated as the source of the dark matter, but as of yet none of these particles have 
been observed. In this paper I present an alternative view. The observations indicate 
the presence of a new force that is much weaker than gravity on small scales and stronger 
than gravity on large scales. This new force is described by a Liouville field that consists 
of a scalar field with an exponential potential. The solutions have a Newtonian behavior 
on small scales and produce forces that are weak compared to gravity. On large scales 
the field has a logarithmic behavior and will therefore produce flat rotation curves. The 
Liouville field may therefore be able to account for the observations that are usually 
taken to imply the existence of dark matter. 
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Introduction 



One of the interesting unsolved problems in modern cosmology concerns the nature of the 
dark matter in our universe. Observations of spiral galaxies show that at large distances 
from the galactic center rotation curves level off instead of having a Keplerian fall off 
|T|. 0, [|. If this behavior is attributed to dark matter we must conclude that about 90 
percent of the mass of spiral galaxies is dark. Observations of the velocity distribution in 
galactic clusters yields similar results. More recently, observations of the CMBR ||, ^], || 
and of type la supernovae |7|, ||, |[ have shown that Qtotai — 1 with Qdm — 0.3 and 
Q DE ~ 0.7, where floM is the density parameter for dark matter and Qde is the density 
parameter for dark energy (also known as quintessence). For excellent reviews of these 



observations see |L0], [LI], [Lg] . 

Of course, all of the above conclusions assume that Newtonian gravity (or general 
relativity) is valid. Another possibility pursued by Milgrom Jl3], |14| and others involves 
modifying Newtonian dynamics. In this theory a = GM/r 2 is replaced by 

fx(a/a )a = GM/r 2 (1) 

where a® ~ 10 -10 m/s 2 , fi(x) ~ 1 when x » 1 and /i(x) ~ x when x « 1. This allows 
one to reproduce the rotation curves of spiral galaxies and some of the other observations 
attributed to dark matter. 

In this paper I will consider a third possibility: there is a new field that produces 
forces which are much weaker than gravity on small scales and which are stronger than 
gravity on large scales. The specific field that I will propose is a scalar field with an 
exponential potential, that is a Liouville field theory ||15|| . I will show that near a point 
source the field behaves as = (fio+a/r and at large distance it behaves as 4> = — ln(r/r ), 
where a, 0o and tq are constants. The Liouville force is given by ac 2 V0, so that at large 
distances the rotation curves will be flat. At small distances the Liouville force will be 
much weaker than gravity if \aa\ << GM/c 2 . Recently there has been another approach 
fl7] , 18, 16] to the dark matter problem that uses a scalar field and an exponential or 
hyperbolic potential. This approach differs from the one presented here in that their 
scalar field does not directly exert forces on particles. This scalar field is a form of dark 
matter as it modifies the spacetime through its energy-momentum tensor and this in 
turn effects the motion of particles in the spacetime. 

Exponential potentials appear in non critical string theory and in dimensionally 
reduced spacetimes with a cosmo logical constant. It will be shown in the last section 
that these theories do not produce Liouville fields with the appropriate parameters to 
be viable candidates for the theory investigated here. 

Throughout this paper I will use the signature (-+++). 
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Field and Particle Equations of Motion 

Consider a collection of massive particles interacting with a Liouville field 0. The action 
will be taken to be 



S = -S„m n / J-g^mU-dr n + iS n / A n (r n ) [g^K + Ife 

(2) 

+aE n m n / 0K(r n )]rfr n - ^ /[V^V M - A 2 e 2 *] ^d 4 * 

where x^(r n ) and U^(r n ) are the position and velocity of the n th particle, r n is the proper 
time along its worldline, m n is its rest mass, A„(r n ) are Lagrange multipliers and a, /3, A 
are constants. Note that and a are dimensionless and f3 has the dimensions of m/kg. 
The field equations are obtained by varying the action with respect to <p(x) and are 
given by 

□0 + A 2 e^ = -^S n m n / <S V - x^r n ))dr n (3) 
v ^ 

where □ = V M V M . The equations of motion for the particles are found by varying the 
action with respect to x%(r n ) and A n (r n ). They are given by 



Contracting with [/" gives 



dX 



+ = am n Vy. (4) 

dr n 



- = -am n — . (5) 



Thus, A = —am n (p and the equations of motion for the particles are 

4- [(i - <*4>)K)\ + (i - afir^uzu* = «V>. (6) 

In the Newtonian limit (and assuming that \a<p\ « 1) the equations of motion become 

^ = -V0jv + ac 2 V0 (7) 
where 0at is the Newtonian potential and I have explicitly included factors of c. 

Solutions of the Liouville Field Equations 

In this section I will examine the spherically symmetric solutions of 

V 2 + k 2 e 2<t> = -a(3M5 3 {r). (8) 
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There are two solutions of interest. If <fi is sufficiently negative the e 2 * term will be 
negligible and the approximate solution will be 

= 00 + - (9) 

r 

where a and 0o are constants. Integrating equation (|8[) over a small spherical volume 
gives a = a/3M/4n. For this to be the solution near the origin (f> << and this implies 
that a/3 < 0. However, if we do not allow true point particles in the theory the vacuum 
solution only extends down to the radius of the object, not to r — > 0. Thus, if 0o is 
sufficiently negative then (f> « outside the object even if af3 > 0. So for a/3 > a 
particle will refer to an object whose size is much smaller than the dimensions involved 
in the problem under consideration. Note that for extended objects there will be a 
correction to a = a/3M/47r produced by the e 2 ^ term. For simplicity I will assume that 
this correction is small. 

The second (exact) solution is 

= -ln(Ar). (10) 

This solution is source free even though it diverges at the origin. To see this integrate 
(§) over a spherical volume of radius R, 

lim f V0 • ds + A 2 lim f e 2lp dV = -a(3M. (11) 

R — >0 J R — >0 J 

Since the left hand side vanishes it follows that M=0. The force on a test particle of 
mass m corresponding to this solution is 



amc 2 



r. (12) 



For this to be attractive we require that a > 0. Thus, if af3 > we have (3 > and if 
a(3 < we have /3 < 0. 

These two solutions, @ and (0), are important because numerical integration of 
(§) shows that the solutions are of the form (||) near the origin and are given by (10) at 
large r. This behavior is easy to show for the logarithmic solution. Let 



The perturbation e(r) satisfies 



ln(Ar) + e(r). (13) 



9 d 2 e de , . 

r V +2r * + 2e = (14 > 



to lowest order in e. The general solution is a linear combination of 

-1/2 

' i 



r ' cos I —y- ln(r) I (15) 
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and 

r- 1/2 sin (^-ln(r)\ . (16) 

Thus, if we start near the logarithmic solution and integrate to large r we will approach 
the solution. If we integrate to small r we will move away from the solution. 

Now consider the observational constraints on the parameters in the Liouville theory. 
Since the rotation curves of spiral galaxies are approximately flat for at large distances 
the Liouville field must be logarithmic at these distances and the force it produces must 
be stronger than the gravitational force. Equating the centripetal and Liouville forces 
gives 

v 2 

a = — ~ 3 x 1(T 7 , (17) 

cr 

for v = 150km/s. It is important to note that, for the spherically symmetric solution, 
the logarithmic Liouville solution is independent of the central mass. This implies that 
the rotation speed in the flat region should be the same for all galaxies. Most galaxies do 
have asymptotic rotation speeds between 100-300 km/s, but they are not all the same. 
The above model is however a little oversimplified. Spiral galaxies are not spherically 
symmetric and this will induce a mass dependence in the asymptotic speed. In fact, for 
the cylindrically symmetric case (see appendix) 



v 2 = a 



2tt 



(18) 



where A is the mass per unit length (or the mass in two dimensions). This mass de- 
pendence will however decrease with distance. Big bang nucleosynthesis calculations 
indicate that there is a significant amount of baryonic dark matter. There is not enough 
of it to account for galactic rotation curves (Q ga i ~ 0.1 and Qbar ~ 0.04), but it may 
effect the asymptotic value of the rotation speed. Another possible explanation is that 
the potential needs to be modified. 
In the region where 

af3M . . 

= 0o + ^r (19) 

the force on a particle of mass m is given by 

We therefore require that a 2 \/3\c 2 << 4nG. This gives 

\0\ « 10~ 13 m/kg. (21) 

Below equation (0) I noted that there will be a correction to a = aftM/An for extended 
objects. If this correction is not small for macroscopic objects the above constraint will 
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be changed. This correction depends on the interior solution and is therefore difficult to 
compute. For simplicity I will assume that it is small and adopt the above constraint. 

A problem arises if the field has a logarithmic behavior in the solar system. To 
see why this is the case note that the radius at which the Liouville force equals the 
gravitational force is given by 

r = ^ ~ 4 x 10 9 to. (22) 
ac 2 

This is obviously in conflict with observations. Thus, near the edge of spiral galaxies the 
Liouville field must exhibit logarithmic behavior but within our solar system it must not. 
I will define three regions, a region in which the field behaves as <fi oc 1/r, a transition 
region, and a region in which <j) oc ln(r). 

The distance to the transition region and to the ln(r) region depends on a and 0o in 
the expression <fi = + a/r. As a simple example consider the case a = and 0o << 0. 
It is easy to show that the distance to the transition region is given by 

r ~ -e"* . (23) 
A 

Numerical analysis shows that the distance to the ln(r) region is given by 

r ~ ^e"* . (24) 

For a point particle a = af3M/An and O is undetermined without some additional 
information such as the value of the field at some point. As discussed earlier there 
will be corrections to a = aftM/An for extended objects. However, just as in the 
point particle case the constants a and 0o cannot, in general, be determined without 
some additional information. This occurs since it is not possible to impose a boundary 
condition at infinity. Thus, without some additional information it is not possible to 
predict the distance to the transition region in our solar system and galaxy. 

The Energy-Momentum Tensor 

The energy-momentum tensor for the 4> field, 

T = 4|^> ( 25 ) 



is given by 

If << we have 



rpfjbV 



(26) 
(27) 
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which is the energy-momentum tensor for a Klein Gordon field if (3 > 0. If (3 < the 
field violates the weak energy condition. For = O + oifiM/Anr the energy density is 

a 2 (3M 2 

P 32vr 2 r 4 ' V ; 

which is much smaller than the Newtonian gravitational energy density. 

In the region where = — ln(Ar) it is easy to see that V Q 0V a — A 2 e 2< ^ = and 

T, lu = -V^cjN^. (29) 



Thus, p = and P r = T rr is given by 



The spacetime will not be asymptotically flat due to the slow fall off rate of the pressure. 

It is interesting to note that violations of the weak energy condition can actually 
occur even if (3 > 0. To see this note that the total energy- momentum tensor is given 
byH3 



= S n ^| J {l- a <P)U^ n 5\x-x n {r n ))dr n +^ [v^V0 - (V Q 0V a - A 2 e 2 ^ 

(31) 

For the solution <fi = — ln(Ar) 



^00 •• m 



S n -^(l + aln(Ar))(C/°) 5\x - x n (r n ))dr n . (32) 
V 9 



Since a > 0, T uu will violate the weak energy condition if particles are placed at Ar < 

Throughout this paper I have assumed that the spacetime curvature produced by 
is negligible. If this is not the case the Liouville field would act as dark matter and 
these effects would have to be taken into account. Here I am taking the point of view 
that the Liouville field is not dark matter, but instead directly exerts forces on objects. 
Since the energy density in the = (po + a/r region is much smaller than the Newtonian 
gravitational energy density I will focus on the = — ln(Ar) region. In this region p = 
and P r will be the source of the gravitational field. In the weak field limit 



□ V = -16ttG (T m „ - - Vfxu T) (33) 

where = r/^ + h^ u and T = rj a pT al3 . Let R be the radius at which the logarithmic 
behavior begins for the Liouville field, and let M be the mass of the galaxy. The Liouville 
field will have a negligible contribution inside R (I will neglect the transition region) and 
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the visible matter in the galaxy will have a small contribution outside R. Solving (p3|) 
with h 00 = 2GM/c 2 r just below r = R gives 

2GM 8nG , fr\ 8nG / R\ 

Thus, the gravitational force on a particle of mass m at r > R is 

-t AnGm „ 



(35) 



(3r 

and for the logarithmic Liouville force to dominate we require that 

1/31 » ^r- ~ KT 20 m/% (36) 

This inequality can be relaxed if we allow the Liouville field to act as a source of dark 
matter. Note that |/ioo| = 1 at r ~ Rexp[\/3\c 2 /8ttG], so that the weak field limit is 
valid out to large distances. 

Thus, the final constraints on a and /3 are 

a ~ 3 x 1(T 7 (37) 

and 

KT 20 « |/3| « l(r 13 m/£;c/ (38) 

with the lower limit on (3 only being necessary if the Liouville field is not to act as a 
source of dark matter. The upper limit is based ont he assumption that a = aftM/An, 
which is valid for a point particle. If the corrections to a are significant for macroscopic 
objects this limit will change. 

Liouville potentials occur in non critical string theory non critical string theory. In 
the Einstein frame the effective action for the dilaton in string theory is [EH] 



V ^-(°- 2 )P6- g W~2) 



(39) 



where k 2 = 8ttG and D is the dimension of the spacetime. The action is of the form (0) 
for 2 < D < 26 with 



K 2 



P = 7 7 , (40) 



, 2(26 — D) 



and 



It is eay see that (3 does not satisfy (|3"8"D, so that the dilaton would be a form of dark 
matter. Secondly, (ct') 1//2 is of order of the Planck length, whereas A -1 is expected to 



A = r , (41) 

3a (D - 2) ' K J 
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be much larger. Thus, the string dilaton is not a likely candidate for the Liouville field 
discussed in this paper. 



Liouville potentials can also occur in dimensionally reduced spacetimes ||21|| . In this 
scenario A 2 is related to the cosmological constant so that it could be astronomical in 



scale. However, (3 violates the lower inequality in ( p8|) so that it would act as a source 
of dark matter as well as exerting a direct force on objects. 



Conclusion 

In this paper I examined an alternative explanation to the observations that are usually 
taken to imply the existence of dark matter. In this approach there is a new force 
produced by a Liouville field. On small scales the force is weaker than gravity but on 
large scales it is stronger than gravity and produces flat rotation curves. This Liouville 
field may therefore eliminate the need for a significant amount of dark matter in the 
universe. 



Appendix 

In this appendix I will briefly examine some of the relevant aspects of the two dimen- 
sional circularly symmetric (or cylindrically symmetric in three dimensions) solutions to 
Liouville's equation. The advantage of working in two dimensions is that exact solutions 
are known. A two parameter solution to 



is given by |22 



-In 



V 2 + A 2 e 2 ^ = 
A 



2n 



(42) 



(43) 



where a and n > are arbitrary constants. Including a point source (replace the right 
hand side of ([42]) by — aj3M8 2 (r)) fixes n to be 



n 



1 - 



a(3M 
2vr 



(44) 



if a(3M/2-K < 1. The constant a is left undetermined. In three dimensions the mass M 
becomes the mass per unit length. 

For r « a the acceleration is given by 



a(n — l)c 2 



and for r >> a 



ac 



~apM 
2vr 



a 2 pMc 2 
2nr 

f 

- 2 - . 

r 



(45) 
(46) 
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Note that for sufficiently large or small r the acceleration is independent of a. However, 
for r ~ a the acceleration will depend on a. For gravity to dominate at small r we must 
have a 2 \(3\c 2 « 2uG. 

Finally, I want to consider how a circular shell effects the behavior of a system located 
its center. If a source of mass M at the origin of the shell the interior solution is 



J int 



In 



A 

— i 
2n 



r_y fax 
aiJ \r 



with 



n 



apM 
2vr 



(47) 



(48) 



Note that this is the same expression for n that one would obtain without the shell. Thus, 
the shell can only change the value of the constant a. Since the force is independent of 
a for small enough r, the shell will not effect the motion of objects close to the source. 
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